COMPUTER ANIMATION OF CONIUGATE SURFACES 


by 

P. S. AVADHANI 


CS£ 

iS92 

n 

Vfl 

com 


HI 

C / I 0) 

% 



DEPARTMENT OF COMPUTER SCIENCE & ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANFL 

SEPTEMBER 1993 


COMPUTER ANIMATION OF CONJUGATE SURFACES 


A Thesis Submitted 

In Partial Fulfillment of the Requirements 
for the Degree of 
MASTER OF TECHNOLOGY 


by 

P.S.Avadhani 


, . V ’ * ’ * ? f I 

to the 

Department of Computer Science and Engineering 
INDIAN INSTITUTE OF TECHNOLOGY KANPUR 
September, 1993 




AVA - 


CCB' 


0 


•, ! " i-im'7 

.■\»‘JH 


4m. Na. A. | 


') '? ■./ , j/! 

r** 



CERTIFICATE 




I - I j.i 



It is certified that the work contained in the thesis entitled COMPUTER ANIMA- 
TION OF CONJUGATE SURFACES by P.S.Avadhani, has been carried dut under 
my supervision and that this work has not been submitted elsewhere for a degree. 



Dr.Sanjay G.Dhande, 

Professor, 

Department of Computer Science and Engineering, 
LIT. KANPUR. 


September, 1993 



Dedicated To 
Dura 



ACKNOWLEDGEMENTS 


I would like to express my deep gratitude to my guide, Dr.Sanjay G.Dhande for 
his expert guidance and constant encouragement throughout this thesis work. 

I am grateful to Mr.K.P.Karunakaran, who, despite his busy schedule, always 
helped me whenever I was in trouble, both in academic and non academic matters. 

I would like to place on record my sincere gratitude for the authorities of Andhra 
University, Visakhapatnam for sponsoring me to this program without which this 
work would not have been possible. 

I am grateful to Prof. S. Biswas, Coordinator, CMP for his help. 

I would like to thank the lab staff of CSE department and the software engi- 
neers at the computer center for their cooperation during this work. Special thanks 
to Mr.N.V.Brahmaji Rao, Software Engineer, CSE department, for his help in the 
preparation of this document. 

I thank my friend Mr.M.V.R.Gopal for his helping me in debugging my programs. 

I am thankful to my friends, Mr.P.V.M.Rao, J.venkatesh, Gagan Saxena and 
B.S.S.Pradhan for the fruitful discussions I had with them. 

I am grateful to Mrs. S. Biswas of CAD-P and to my friends Mr.P.V.S.N.Murthy, 
Mr.S.P.Chandrasekhar, Dr.V.C.V.Rao and Mr.Appaji Rao for their company. I am 
thankful the staff of the CAD-P for their help . 

My special thanks to Mr.P.V.Ravindra and Mrs.Sunitha Ravindra for the delicious 
Sunday dinners I had in their house during the last two months of my stay . 

I could not have completed this work without the support of my wife, Lakshmi, 
who, despite her own domestic troubles kept on encouraging me. 


— P.S.Avadhani. 



Contents 


1 INTRODUCTION 1 

1.1 CONJUGATE SURFACE GEOMETRIES 2 

1.2 COMPUTER ANIMATION 5 

1.3 LITERATURE REVIEW 8 

1.4 SCOPE OF THE PRESENT WORK 10 

1.5 ORGANIZATION OF THE THESIS 11 

2 COMPUTATIONAL CONJUGATE GEOMETRY 12 

2.1 SYMBOLIC MODEL 12 

2.2 BASICS AND OVER VIEW OF THE SYMBOLIC ALGORITHM . 12 

2.3 THE SYMBOLIC ALGORITHM 18 

2.4 EXAMPLES 20 

2.4.1 TWO DIMENSIONAL CAM 21 

2.4.2 CYLINDRICAL CAM - TRANSLATING FOLLOWER .... 24 

2.5 ANIMATION AND GRAPHICS DISPLAY 25 

2.6 ALGORITHM 26 

3 COMPUTER ANIMATION OF CONJUGATE SURFACES 30 

3.1 INTRODUCTION 30 

3.2 BASIC CONCEPTS 30 

3.2.1 DEFINITION OF A GENERIC CUTTER 31 

3.2.2 DESCRIPTION OF THE GENERIC CUTTER 31 

3.2.3 ALGORITHM 35 

3.3 METHODOLOGY 37 


V 



vi 

3.3.1 ALGORITHM 38 

3.3.2 ALGORITHM 39 

3.3.3 ALGORITHM 41 

4 SOFTWARE DEVELOPMENT 44 

4.1 CUTTER PROFILE MODULE 44 

4.2 SYMBOLIC COMPUTATION MODULE 44 

4.3 SURFACE GENERATION MODULE 45 

4.4 CUTTER SURFACE GENERATION MODULE 46 

4.5 TRIMMING SURFACE GENERATION MODULE 48 

4.6 ANIMATION MODULE 48 

5 CONCLUSIONS 49 

5.1 SCOPE FOR FUTURE WORK 49 

APPENDICES 51 

A MACSYMA 51 

A-1 CAPABILITIES 51 

B STARBASE GRAPHICS LIBRARY 53 

B-1 DOUBLE BUFFERING 53 

B-2 TRANSFORMATION MATRIX STACK 53 

B-3 HIDDEN SURFACE REMOVAL 54 

B-4 RENDERING 54 

C BIBILIOGRAPHY 56 



List of Figures 

1-1 Definition, Display and Realization of a Circular Shape 3 

1- 2 Definition, Display and Realization of a Straight line Shape 4 

2- 1 A Pair of Conjugate Surfaces 14 

2-2 Translational and Rotatioinal Degrees of Freedom of Conjugate Coor- 
dinate Systems 15 

2-3 Two Dimensional Cam - Translating Flat faced Follower 23 

2-4 Cylindrical Cam - Translating Follower 23 

2-5 Animation of Two Dimensional Cam - Flat faced Follower 28 

2-6 Conjugate Surface - Three Dimensional Cam 29 

2-7 Animation of Conjugate Surface - Three Dimensional Cam 29 

2- 8 Animation of Conjugate Surface - Three Dimensional Cam and Follower 29 

3- 1 Definition of a Generic Cutter 34 

3-2 A Particular Cutter Surface 42 

3-3 Initial and Final Trimming Surfaces 42 

3-4 Animation of a Conjugate Surface ( Translating 50 units along X-axis ) 43 

3- 5 Animation of a Conjugate Surface ( Translating 50 units along X-axis 

and Rotating 15° around X-axis) 43 

4- 1 Software Development 47 



ABSTRACT 


Computer Animation has been a topic of research for many scientists involved with 
human body modelling and movement, education and visualization of mathematical 
and scientific models. The main thrust of most of the research has been towards the 
human body modelling. 

The report proposes a methodology for animating two moving surfaces in higher 
pair contact (called the conjugate surfaces). The methodology is ba.sed on a sym- 
bolic algorithm which has been implemented using a symbolic computation program 
MACSYMA. Two case studies, namely, a 2D cam and a 3D cam with translating 
followers have been shown in animation using this algorithm to illustrate the model. 

A general model of the cutter is proposed and the surface generated by the cutter 
is shown in animation from a specified initial point to a specified final point. Linear 
interpolation is used as the in-betweening technique. 

The present work can be applied in Design and Manufacturing field as it illustrates 
the design of the cutter surfaces. 

This work is implemented using the starbase graphics library on HP-9000/834 
Turbo SRX 3D graphics work stations. 



Chapter 1 


INTRODUCTION 


Recent developments in Computer Graphics have opened new dimensions in the fields 
of Animation, Graphics, Design and Scientific Visualization. Primarily for animation, 
handling of huge data sets and the speed are the basic problems. This is to some ex- 
tent being solved at the advent of the high speed computers with high data handling 
capabilities. It is also to be noted here that the geometry of the object or scene to 
be displayed in animation plays another important role. The motion dynamics[16], 
the in-betweening techniques[22] add to the realism in animation [43]. For example, 
surface generation is a major area of interest in solid modeling, design and manu- 
facturing science. Suppose, a cutter generates a surface while cutting a blank, then 
to show that surface with cutter in animation, we should know the geometry of the 
cutter, the geometry of the surface to be generated and the mode of in-betweening 
apart from the relative motion specifications. In fact, the developments in geometric 
and solid modeling techniques have provided the required base for these problems 
for design and manufacturing engineers with an elegant and precise definition of a 
component or product that needs to be designed[31,64,67]. In manufacturing science 
several methodologies based on these modeling techniques were developed. Recently 
a unified methodology for geometric modeling of various surfaces using the concept of 
Conjugate Geometry is developed. This methodology is useful in developing any two 
surfaces in higher pair contact and in movement. The work proposed in this thesis 
is to display the cutter from a specified initial position to a specified final position 
generating the surface in motion using the symbolic algorithm based on Conjugate 
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Geometry. The following sections illustrate the basic concepts of Conjugate Geometry 
and animation. 

1.1 CONJUGATE SURFACE GEOMETRIES: 

It can be seen from the experience that a mathematical description of the shape 
of an object will be specified depending on the purpose of its use. For example, if a 
shape is to be displayed on a computer controlled display device then a parametric 
definition is appropriate whereas if a shape is to be specified for some analytical 
procedure then a closed form expression in Cartesian coordinates is used. However, 
if a shape is to be realized by generating process then neither of these is useful. 

Consider the case of generating a circle from a 2D lamina. One way could be to 
employ a series of straight line cuts (see Fig.1-1). These cuts constitute a family of 
curves. An envelope to this family is the shape generated , which in this case, is a 
circle. Similarly, a straight line can be realized as an envelope of a family of circles. 
The family is generated by moving the center of the circle on a linear path. Thus, 
a family of circles generates a linear shape and a family of lines generates circular 
shape ( See Fig. 1-2). Thus envelope theory is the mathematical basis for modeling 
the geometry of several shape generating processes. Traditionally the envelope theory 
defines an enveloping curve or surface in a fixed frame of reference. However, many 
manufacturing processes require moving frames of references. The conjugate geometry 
model [28] extends the classical envelope theory and defines the envelope in one 
moving frame of reference due to a family generated by another shape ( curve or 
surface ) and the relative motion between the two frames (see Fig.2-1). The present 
work uses the conjugate geometry approach to simulate the animation of a cutter to 
generate a surface from a given initial position to a given final position. 




Figure 1-] Definition, Display and Realization of a Circular Shape 

(a) Mathematical Definition - Implicit Form 

(b) Display Definition - Parametric Form and 

(c) Realization Definition - Envelope Form 






Y 


(C) 



Envelope of 
a family of circles 


Figure 1-2 Definition, Display and Realization of a Straight Line Shape 

(a) Mathematical Definition - Implicit Form 

(b) Display Definition - Parametric Form and 

(c) Realization Definition - Envelope Form 
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1.2 COMPUTER ANIMATION: 

Animation means bringing to life. In this view, people often think of animation 
cLS synonymous with motion. However, it covers all changes that have a visual effect. 
Thus, it includes the motion dynamics ( i.e. the time varying positions ) and update 
dynamics (shape, color, transparency and texture ) and all changes of rendering effects 
of an object. Animation is used widely in the entertainment industry. At the advent 
of computers, it is also being applied in education, in industrial applications such as 
control systems, flight simulators for aircraft and in scientific research. The appli- 
cation of Computer Graphics and especially of animation have come to be grouped 
under scientific visualization. Often, the animations in scientific visualization are 
generated from simulations of scientific phenomenon. The results of the simulations 
may be large data sets which are converted into images that constitute the animation. 
Visualization is more than mere application of graphics to science and engineering 
and it can involve other disciplines such as Computational Geometry and Geomet- 
ric modeling. The present work is to show a manufacturing process in animation. 
All manufacturing processes are shape realization processes. The shape realization 
is equally important as its design. The proposed animation enables a manufactur- 
ing engineer to evaluate the likely outcome of a particular cutter of a manufacturing 
process. 

Conventional animation is oriented towards the production of two dimensional 
cartoon films and it is a fairly fixed sequence. The steps involved in this are 

1. Story: 

a ) Synopsis - The summary of the story in a few lines. 

b ) Scenario - Detailed text that describes the complete story without any cine- 
matographic references. 

c ) Story board - This is a film in outline form. It consists of a number of 
illustrations arranged in comic-strip fstshion with appropriate captions. The 
number of individual illustrations within a story board will vary widely. The 
important aspect here is that they represent the key moments of the film. H 
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is also important to note that the film is composed of sequences that define 
specific actions. Each sequence consists of a series of scenes that are generally 
defined by a ertain location and set of characters. Scenes are divided into shots 
that are considered as picture units. 

2. Layout - This step consists of the design of the characters to be animated and 
action plotting. Based on the story board the relationship between the shapes 
and the forms in the background is decided. The layout department in an 
animation studio has to finish drawings of the settings and sketch background 
layouts. The artists must have a knowledge of the physical characteristics of 
the camera that will be used to photograph the frames. 

3. Sound track - In conventional animation, sound track recording has to precede 
the animation process because the motion must match the dialogue and /or 
music. 

4. Animation - The animation process is carried out by the animators who draw 
key frames. Often aji animator is responsible for one specific character. 

5. In-betweening - ’’In-betweens” are defined as drawings which are placed between 
two key positions or frames. Assistant animators draw some in-betweens and 
in-betweeners draw the remaining figures. The work of the assistant animators 
requires more artistry than the in-betweeners, whose task is almost automatic. 

6. Xeroxing and Inking - Sketches are usually drawn and then they have to be 
transferred to acetate cels using modified xerox cameras. Lines will be inked in 
by hand. 

7. Painting - As cartoon animated films are usually in color, they must go through 
a painting stage. This work requires patience and accuracy. Cels must have the 
right degree of opacity and static backgrounds must also be painted. 

8. Checking - Animators need to check the action in their scenes before shooting. 
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9. Shooting - The final photography of composite animation is usually done on 
color films or video tapes. The shooting phase is not a trivial operation. Move- 
ment can be simulated at this stage by moving certain cels with respect to 
some others. Many cels are placed in layers one above the other and are shot 
simultaneously. 

10. Editing - This step is considered part of the post production stage. 

Introduction of computer in the field of animation came as more than automating 
the various steps in the traditional animation process. The initial application was in 
the in-betweening, coloring and copying phases. But soon it developed into an entirely 
new field by itself. This was brought about mainly by the ability of the computers 
to represent, process and operate on objects of dimensions higher than two. Now the 
main task lies in modeling, representing and simulating dynamic phenomena of all 
kinds. 

Most of the phenomena which may be represented on the screen of a workstation 
are typically time-dependent. The techniques of Computer Graphics allow construc- 
tion of 3D graphical objects using geometric modeling techniques. Moreover, in a 
3D space, scenes are viewed using vertical cameras and they may be lit by synthetic 
light sources. In order to visualize these phenomena at any given time, it is neces- 
sary to know the appearance of the scene at this time and then computer graphic 
techniques allow us to build and display the scene according to viewing and lighting 
parameters. The general problem in this situation is to express time dependence in 
the scene and how to make it evolve over time. These problems and their various 
solutions are part of what is usually understood by the term Computer animation. 
Thus computer animation may be defined as a technique in which the illustrations of 
movement is created by display on a screen, a series of individual states of a dynamic 
scene. Broadly, there are three ways of animating this ,namely, rotoscopy, key frame 
animation and procedural animation. 

In key frame animation, a certain number of frames called key frames are given to 
the system and the computer derives the other frames using interpolation procedures. 
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Various approaches have been developed in this process. The two most important and 
widely used techniques are shape interpolation and parametric key frame animation. 
In shape interpolation the in-betweens are obtained by interpolating the key frame 
images themselves. In parametric key frame animation objects are characterized by 
parameters; motion is specified by giving key values for each parameter; in-between 
values are calculated using an interpolation law. In procedural animation, motion is 
algorithmically described. 

Initially and even now, most of the in-betweening work is done using the technique 
of linear interpolation. However, it is observed that this might cause a lack of smooth- 
ness which considerably alters the motion. Undesirable effects such as discontinuities 
in the speed of motion and distortions in rotations are produced. Consequently for 
such motions, other methods should be used. A non-linear interpolations (non linear 
in time but not in space ) like spline interpolations are some of the methods tried by 
researchers with good results. 

1.3 LITERATURE REVIEW: 

The idea of animation is conceived as early as 1831 by Joseph Antoine Plateau 
who produced a device called phenakisto scope. This consisted of a spinning disk 
with drawings that could be viewed through framing windows. Almost at the same 
time, Horner produced a device called the Zoetrope, in which drawings were produced 
on the inner surface of a drum. The individual frames could be viewed through slits 
in the inner rotating drum as they came past a fixed framing slit, through which the 
viewer looked. Devices of this sort became highly popular during the last half of the 
previous century. 

Paralleling this development in frame animation, photography evolved fast in 
the last half of nineteenth century giving boost to the field of animation. In 1915, 
Earl Hurd introduced the technique of cel animation which took its name from the 
transparent sheets of celluloid that is used. But, of course, the father of commercial 
animation is certainly Walt Disney . In the ten years 1928-1938 he produced Mickey 
Mouse Donald Duck and the Silly Symphony series. 

The animation film industry was well established when the computer made its first 
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contributions. One of the early pioneers were Bell Laboratories. This work was pri- 
marily directed at the scientific and education world. However, it was not long before 
the importance of the computers to the animation industry was realized. Advances 
in real-time systems allowed the computer to automate the control of the rostrum 
camera. Kallis[38] in 1971 described a system capable of carrying out sequences of 
operations using pan, tilt, spin and zoom. Following this lot of work has been done 
in the automation of in-betweening. 

Reeves [79] proposed a method called moving point constraint method. This 
method is meant to allow the specification of multiple paths and speeds of interpola- 
tion and to reduce motion discontinuities at key frames. The principle of the technique 
is to associate a curve varying in space and time with points of the animated objects 
called the moving point and it controls the trajectory and the dynamics of the point. 
Kochanek-Bartels[51] used spline interpolation to control the tension, continuity and 
bias and thereby controlling the in-betweens. These methods coupled with methods 
based on aerodynamics are used to simulate and control the motion in objects in fluid 
flows [95]. Freedman et al[35] used a technique of quadrature pair of oriented filters to 
vary the local phase, giving the sensation of movement for displaying patterns that 
appear to move continuously without changing their positions. 

An interactive technique for animating deformable objects is developed in [26]. 
This approach provides a representation of the deformations independent of the sur- 
face geometry and can be easily integrated into traditional hierarchical animation 
systems. 

A general interactive modeling and animation system that facilitates a variety of 
simulation and animation paradigms like cameras, light sources, rendering and user 
interfaces have been proposed in [101]. Much work has been done in human body 
modeling and the pioneers in this area are Thalmann [59,60,61,62,63]. 

Human hand motion and grasping of arbitrary shaped objects is developed in 
[82] using a knowledge based approach. This effectively reduces the enormous search 
space associated with the problem of grasping objects and satisfying the kinematic and 
physical constraints associated with the problem. A methodology for controlling the 
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dynamic legged motion is proposed in [77]. LcLsseter[54] ,Laybourne[56] and Shoup[88] 
give a very good account of the techniques of traditional animation applied to 3D 
Computer animation apart from the various other techniques of in-betweening and 
rendering. Burtnyk[16] proposed techniques to enhance the motion in key frame 
animation. 

Animation of facial expressions was dealt by Parke[74], Platt and Badler[75], Thal- 
mann[58] and Max[65]. Shoemake [87] proposed a method based on quaternions for 
rotations in animation instead of conventional methods. 

The developments in geometric and solid modeling techniques [67, 64] have given 
the manufacturing engineer the required methodology. The problem of conjugate 
geometry between a pair of moving as well as contacting surfaces is termed as a higher 
pair contact problem in kinematics[7, 23, 97]. Such problems have been analyzed by 
designers of cams and gears. Chakraborty and Dhande[23] described the methodology 
of defining the kinematics of relative motion as well as the geometry of the contacting 
surfaces for three dimensional cams. One of the major computational difficulties in 
formulating and analyzing the problems of conjugate geometry is that for every given 
situation of a pair of contacting surfaces and their relative motion, one has to derive 
the biparametric equations of the surfaces every time. This can be improved by using 
the symbolic manipulation approach [28, 45]. A symbolic model with twelve degrees 
of freedom is proposed in [28]. 

1.4 SCOPE OF THE PRESENT WORK: 

The idea of the present work is to depict the power of conjugate geometry approach 
to display conjugate surfaces in animation. Two examples of 2D and 3D cams with 
translatory motion are displayed in animation using the symbolic algorithm. The 
next phase can be viewed as two parts: The surface generation and the animation. 
A generic cutter with nine parameters is defined and displayed. The cutter draws 
the surface it generates from a given initial position to a given final position while 
moving along the path. The initial and final positions will have three rotations and 
three translations as their input. It means that the cutter can translate and/or rotate 
along any axis while it is moving. 
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Initially, the cutter is drawn with the given nine paranrieters( which determine 
the cutter) using an algorithm. Next, the symbolic program MACSYMA is used to 
generate the C-code for the conjugate surface generated by the cutter. In fact, it is a 
set of eight surfaces that are drawn, as the cutter shape has eight segments, namely, 
six straight line segments and two circular segments. Next, a trimming algorithm 
is used for depicting the ends of the surfaces generated. Thus the total number of 
surfaces generated will be ten in number. These ten surfaces are properly combined 
along with the cutter itself and are generated part by part for each time frame to give 
the effect of continuous motion. 

All these algorithms are implemented using the double buffering technique of 
starbase graphics library on HP-9000/834 Turbo SRX 3D graphics work stations 
based on HPPA RISC CPU rated at 2.02 MFLOPS(14 MIPS) performance. 

1.5 ORGANIZATION OF THE THESIS: 

In chapter 2, a detailed account of the computational conjugate geometry along 
with the symbolic algorithm are presented. In section 2.4 two case studies, namely, 
the 2D and 3D cams are presented to illustrate the power of the symbolic algorithm. 
The animation aspect is also taken into account. Chapter 3 deals with the generic 
cutter definition, the algorithm to generate the cutter and the methodology cidopted 
to generate the surfaces and animation. Chapter 4 deals with the aspects of software 
developments and methodology adopted and finally chapter 5 offers conclusions. 



Chapter 2 


COMPUTATIONAL CONJUGATE GEOMETRY 


2.1 SYMBOLIC MODEL: 

There are basically three elements involved in any machining process. The first 
one is the surface of the cutter denoted by Ei ,the second being the surface that is 
to be produced denoted S 2 (see Fig.2-1). The third element is the relative motion 
between these two surfaces which will be constrained by the type of kinematic chain 
of the machine tool. The surfaces Si and S 2 are conjugate to each other as those two 
surfaces maintain a higher pair contact throughout the cutting motion . A unified 
symbolic model based on the computational conjugate geometry is proposed in [28] 
with twelve degrees of freedom. Using this symbolic model, if any two of the three 
elements, namely, the surfaces Ei, E 2 and the relative motion parameters are known, 
the third one can be determined. This model is useful in finding the NC cutter path 
when the surfaces Ei and S 2 along with the relative motion parameters are required 
to be determined. This model is useful in NC path simulation, NC code optimization 
etc. The best use of this model may be for choosing the optimal combination of the 
machine tool and the cutting tool by an iterative process. 

2.2 BASICS AND OVER VIEW OF THE SYMBOLIC ALGORITHM: 

A geometric model of a machine tool is essentially a representation of the kinematic 
structure of the machine tool . A kinematic structure of a machine tool describes the 
kinematic chain necessary for providing the cutter motion to the cutter and the feed 
motions to the cutter and /or the blank. The symbolic model considers the cutter 
and the blank to be two free bodies in space and thus their position and orientation 
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can be described by six degrees of freedom each. These twelve degrees of freedom 
are expressed with respect to a fixed frame of reference So{0o — Xo,yo, Zq). The 
coordinate systems 51(01 — Xi^Yi^Zi) and ^ 2(02 — A'’ 2 ,y 2 j-^ 2 ) ^re considered 

attached to the cutter and blank respectively (see Fig.2-2). The initial positions of 
Si and S 2 are denoted by Si{Oi — Xi,Yi,Zi) and S^iO^ — X^iY^., Z^)- 








Let the six tuple (a., 6 ., c,, <?,, = 1,2 denote the translations and rotations 

of the moving coordinate frames Si and S2 respectively. These twelve parameters 
of motion represent the kinematic arrangement of the machine tool so that they 
can be expressed as functions of time t. As the tool surface Ei and the surface S2 
produced by blank are time-invariants, the known surface among the two surfaces can 
be expressed as a biparametric surface in u and v. Since the two Sj and E2 surfaces 
are in higher pair contact, at any point of contact P, 

[“Ail = ’? (“Ail 

and 

= ±'2[oM] 

where 

f - Point P with reference to Si 

2 

? - Point P with reference to S2 

- Unit normal vector of the surface Ei at point P with reference 
to Si 

2 

- - Unit normal vector of the surface S2 with respect to S2 

[j M] - Homogeneous transformation matrix transforming a vector from 
the coordinate system Si to 

[2M] - Homogeneous transformation matrix transforming a vector from 
the coordinate system ^2 to 

It is to be noted that all vectors such as position vectors, normal vectors and ve- 
locity vectors are all expressed as row vectors in the homogeneous coordinate system. 
Hence the fourth element of any position vector is unity and it is zero of the normal 
vectors. 

The surface S2 is the boundary of the swept volume of the cutter surface Ei which 
is nothing but the envelope to the various instances of the cutter surface. In order to 
identify the point on the cutter surface Sj which is cutting the blank and producing 
the surface E2 at any given time, it is necessary to have some additional constraint 
called the condition of contact. It can be said that, in order to ensure a smooth 
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contact or cutting at given time time, the relative velocity vector between the cutter 
and the blank must be orthogonal to the common normal at the point of cutting [23]. 
Mathematically, it can be expressed as 


^ P.12 _ Q 


and 


P.12 _ Q 


where 


- - Unit normal vector of the cutter surface Ej at point P with reference to the 
coordinate system So 

- Unit normal vector of the cutter surface S2 at point P with reference to the 
coordinate system So 


Oy 


- - Relative velocity vector between the tool and blank with reference to the 

coordinate system 5 o 

The condition of contact is a scalar equation with parameters u, v and t. Thus 
the condition of contact is a relation between u and v of any given instance of 
time. For any value within its range , a corresponding value of the parameter 
v can be found. Thus one can locate a set of points or a curve on the surface 
which defines the surface curve for E2 as well. The corresponding points of the 
surface can be obtained using the equation 


>1 [jm] 


Formally the symbolic algorithm can be stated as follows: 
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2.3 THE SYMBOLIC ALGORITHM: 

2 

1. Define the position vector ? of a generic point P on the cutter surface as 
a biparametric surface in u and v such that 

V 

= [xi(w,u) yi{u,v) Zi(u,u) 1] 

item Define the twelve parameters of motion {ai,bi,Ci, 6 i, = 1,2 in 

terms of time t. 

2. Calculate the matrices [jM] and JjA/j as follows: 

;?m] = [;r,] [;/i,] [;r] 

where 


”1 0 0 0 " ’cos<^i 0 — sin<jii O' 

0 D ] _ 0 cos^i sin^i 0 foo 1 _ 0 1 0 0 

H 0 —sin 01 cos 01 0 sin^i 0 cosi^i 0 

.0 0 0 ij L 0 0 0 1. 


cos 01 

szn0i 

0 

O' 


■ 1 

0 

0 

O' 

- sin 01 

cos 01 

0 

0 

■|-J o' 

. 

II 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1. 


.Oi 

^>1 

Cl 

1. 


3 . Calculate the matrices ’^M ] , ^M] and ^ follows: 

[2m] = [;b,] [X] [2K.] [Sr] 

where 



'1 

0 

0 

O' 


'cos 02 

0 

— sin 02 

O' 

Op] 

2iXx — 

0 

cos 02 

sin 02 

0 

fo p ] 

[ 2 -^yl — 

0 

1 

0 

0 

0 

— sin 02 

cos 02 

0 

sin 02 

0 

cos 02 

0 


.0 

0 

0 

1. 


0 

0 

0 

1. 
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cos V>2 sint/>2 0 0 

] _ — sini/’2 005 02 0 0 

■I “ 0 0 10 

. 0 0 0 1 



■1 0 0 O' 

0 10 0 

0 0 10 

<12 62 C2 1 . 




-1 


Matrix [ 2 -^] is the time-<iifferential of matrix [^M] 


4. Derive the position vector 


P( . 


[X2(u,l)) Z/2(u,v) 22 (u,v) l] 


using 



5. Derive the normal vector 




p.’ 


du dv 


where 

g' P 3 ' p 1 

and are the partial derivatives of ? with respect to u and v 
respectively so that 


a'? 

— r 


dzi 

du 

“ 1 du 

du 

du 

d'l 

[■ Bx\ 

§}a. 

^£1 

dv 

1 dv 

dv 

dv 


0 ] 


0 ; 




p.i 


i 

j 

k 

dxi 

9yi 

dzi 

du 

du 

du 

dxt 

dvi 


dv 

dv 

dv 


and hence 
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6. Calculate the unit normal vector 


where || is the norm of 

7. Calculate the unit normal vector in the global coordinate system 

using 


0 


n 


p.i 




8. Calculate the relative velocity vector 


where 





and 


p.j _ 


\M 



9. Reduce symbolically the condition of contact 


10. obtained in step 5 will be a function of u, v and t. Use the condition of 
contact obtained in step 10 to eliminate either of u and v. This will yield 

S2 

2.4 EXAMPLES 

In this section two examples of the cams are discussed along with the implemen- 
tation of the symbolic algorithm. Both of these are implemented using starbase 
graphics[47] using double buffering and are shown in animation. 
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2.4.1 TWO DIMENSIONAL CAM: 

The first example is a two dimensional disc cam with a translating flat-faced 
follower (Fig. 2-3). The flat face is perpendicular to the direction of translatory 
motion of the follower; a is the offset distance ana r is the radius of the circular 
cam [23]. 

The cam rotates with the uniform angular velocity t. The angular rotation 
is the parameter of the cam motion. The follower’s linear displacement k is 
measured from Oq. 

The twelve parameters of motion are 

ai = k sint 02 = 0 
fei = — fc cost bi = r — k cost 

Cj = 0 C 2 = 0 

^ 1=0 ^2—0 

,^ 1=0 <^2=0 

i’l = t ^2 = 0 

and 

r • 1 

. = [r cosu r sinu v 1] 

Proceeding as discussed in the section 2.2 above, 

_ = [fc sint + r cos(u -f t) — r -|- r sm(u d- 1) v 1] 

with the condition of contact 

TT 

» + (= 


Eliminating t, we get 
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X2 = :tk cosu 

j/2 = 0 
Z2 = V 

As the discussion is about a two dimensional cam, putting z = 0, the corre- 
sponding curve is given by 


X2 = ilk cosu 
y2 = 0 

Z2 = 0 


which is the equation of E 2 . 

This is implemented in starbase and the cam in motion in graphic display. A 
particular instance is shown in Fig.2-4. 
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2.4.2 CYLINDRICAL CAM - TRANSLATING FOLLOWER: 

The second example considered is a cylindrical cam with translating offset fol- 
lower as shown in Fig 2-4. Here a is the amount of offset and the parameters 
are follows: 


and 


Oi = —a 

bi=0 
s . 

Cl = — smt 

6r = - 
* 2 

<Ai = 0 


= 0 


02 = 0 
62 = 0 

C 2 = 0 
^2 = 0 
(j>2 = 0 
^2 = t 


1 


P 


= [r cosu 


r sinu v 1] 


Thus 


X 2 =r cosu cost — V sint — a cost 


J/2 


-r cosu sint — v cost -f a sint 


Z 2 = r sinu + - sint 

2 


with the condition of contact 


u = — cost tanu 

2 

Eliminating t from the two equations (2.1) and (2.2) the surface E 2 can be 
obtained. 
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2.5 ANIMATION AND GRAPHICS DISPLAY: 

The graphics display of the first example is fairly straight forward as the surfaces 
El and E 2 are a circular disc and two parallel lines touching the circle. (Fig. 2- 
4). However, the display of the cylindrical cam with translating follower is 
a bit tricky. The symbolic algorithm generates only the surface S 2 which is 
conjugate to the cutter surface. To have the realism in display and animation, 
it is necessary to draw the cylindrical cam and the cutter along with the surface 
generated in motion. This requires finding the intersection curves of the cylinder 
with the surface £ 2 - This results in the following discussion: 

Let R be the radius of the cylinder with it’s angular velocity t. Then the 
equation of the cylinder is 


X = R cost 
y — R sint 
z = d 

where d is the height of the cylinder(see Fig 2-6). 

The surface calculated from the symbolic algorithm, given in equation(2.1) and 
(2.2) is 


X 2 = r cosu cost — V sint — a cost 
j /2 = — r cosu sint — v cost -f a sint 


. 

Z 2 = T sinu -f - sint 

JL 


with 


V — — cost tanu. 

2 

The intersection curve of the cylinder and the surface should be a circle in the 
XY-Plane (see Fig.2-5). Hence, 
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-^2 , „,2 J 1 , ,2 

a: -|- 2/ — arj 4 - 2/2- 

Simplifying and eliminating v, we get 


4 cos^u — 8 a r cos^u + {4 — A B? — cos^t) cos^u -f cos^t = 0. 

This is a fourth degree polynomial equation in cosu for each instance of t and 
hence gives rise to four roots for cosu. However, the roots are meaningful only 
if they lie in the interval [—1,1]. 

Analyzing further mathematically with the constraint that |cosu| < 1, it can be 
seen that only two valid roots for cosu are possible for each instance of t. These 
two roots give rise to two curves which define the boundary of the surface S 2 
as shown in Fig.2-5. 

Once the curves are found then the following algorithm illustrates the display 
and the animation of the cylindrical cam along with the translating follower. 
Some instances of these are shown in Figures 2-6,2-7 and 2-8. 

2.6 ALGORITHM: 

1. Define the cutter surface 

[r cosu r sinu v 1] 0 < v < d 

2. For t = 0 to MAXT { 
a ) Define the cylinder 


X = R cosu 


y = R sinu 
z = fi{t) 
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where fi{t) is one of the two intersection curves nearer to z = 0 
b ) Define the cylinder 


X = R cosu 
y — R sinu 
Z = / 2(0 

where is the intersection curve nearer to z = d. 
c ) Define the cylinder 

X = R cosu 
y = R sinu 

2 = U flit)<V<f2(t) 

d ) Define the circular discs 

X = Ri cosu 
y = Ri sinu 
z = flit) 

and 

X = Ri cosu 
y = R\ sinu 
z = fiit) 

} 














Chapter 3 


COMPUTER ANIMATION OF CONJUGATE SURFACES 


3.1 INTRODUCTION: 

In this chapter, the algorithms for generating the cutter surface, it’s conju- 
gate surfaces and the algorithms used for animation of conjugate surfaces axe 
presented. A model of the cutter determined by nine parameters [27] is used 
to generate the cutter. Another twenty four parameters determine the initial 
and final positions of the cutter. These parameters are the twelve parameters 
of the symbolic algorithm for the initial position of the cutter and the twelve 
parameters for the final position of the cutter. 

The motion of the cutter or the in-betweening is simulated through the tech- 
nique of linear interpolation between the initial and final parameters of the 
cutter. The details are discussed in subsequent sections. 

3.2 BASIC CONCEPTS: 

The symbolic algorithm coupled with the technique of in-betweening using lin- 
ear interpolation is the ba.sic methodology used for producing the animation . 
Before going into the details, it is necessary to describe the definition and the 
algorithm to generate the generic cutter. 
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3 .2. ] DEFINITION OF A GENERIC CUTTER: 

The geometry of a cutter can be modeled as an axisymmetric, biparametric 
surface. Dhande and Karunakaran [27] proposed a model of the generic cutter 
in which the generatrix curve is a planar curve consisting of a set of contigu- 
ous straight line segments and circular fillet segments. The directrix motion is 
the axisymmetric sweep of this curve around an axis in the plane of the gen- 
eratrix curve. The generatrix curve is assumed to be a piecewise continuous 
curve consisting of three straight line segments and a circular fillet (see Fig.3- 
1). This generatrix curve can be defined by means of eight parameters, namely, 
d, r, e,/, 0,6, A and hi as shown in Fig.3-1. The cutter surface can be obtained 
by rotating the generatrix curve along the z-axis. 

3.2.2. DESCRIPTION OF THE GENERIC CUTTER: 

The planar curve of the generic cutter surface consists of three straight line 
segments and a circular arc segment(Fig.3-l). A coordinate system Sc{Oc — 
XciYc, Zc) is considered to be attached to this cutter profile. The planar curve 
of a cutter profile can be represented by means of a parameter u which is the 
distance of any generic point on the profile with respect to the origin Oc- The 
cutter surface Si is obtained by sweeping this generatrix profile along Zc axis. 
The four segments of the generatrix curve can be algebraically expressed as 


where 


OcP : Xc tana 


PQ : {xc - eY + {zc - ff = 

d d 

QS : Xc = Zc tanb — - tana tanb -f - 

z z 

ST : Xc = h tanb — ^ tana tanb + - 

z z 


OcP - First straight line segment 

PQ - The circular segment with center at A 
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QS - Second straight line segment 
ST - Third straight line segment 

R - Point of intersection the first segment OcP and the third segment QS. 
d - Diameter of the cutter mea.sured at R 

r - Radius of the circular segment. This is positive for convex arc and 
negative for concave arc. 

a - Angle of segment OcP measured from Xc axis (— | < a < |). 
b - Angle of segment QS measured from Zc axis (— ^ < a < |) 
e - X coordinate of point A in the coordinate system Sc 
f - Z coordinate of point A in the coordinate system Sc 
h - Non parallel height of the cutter 
hi - Parallel height of the cutter 

is-fillet-radius - Flag set to TRUE if OcP and QS are tangents to the 
circular segment PQ else FALSE. 

u - Parameter which defines the distance of any generic point P on the 
profile measured from Oc 

V - Parameter describing the sweep of the two dimensional profile 

[JM] - Homogeneous transformation matrix describing the sweep of the 

two dimensional profile 

•51,52,53,54 - lengths of the segments OcP, PQ,QS and ST respectively 

- position vector of a generic point P on the two dimensional profile of 
the cutter in Sc coordinate system in terms of u 

- Position vector of a generic point P on the surface of the cutter in Si 
coordinate system in terms of u and v 

The parametric equation of the generatrix curve ?(u) with respect to the co- 
ordinate system Sc{Oc — Xc,Yc,Zc) is defined by 

Cp 

^ (u) = 1] 
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whore 


3 -c(u) 


u COS a. 


0 < u < Si 

c + lr| cog { arc! an ) + (^)}, Si < u < Si + 52 
u - (gj + S2 + 53) sinb + X*, 




5i + ^2 <C ti 

“h ^2 ^3 ^ *^2 ^3 


ydu) - 0 




{u) = { 


u sina^ 


0 < u < $1 


f + |r| gm{arctan(- ^»-^ - _ -f) + (^)), 5 i < u < Si + 52 

It — (5i + 52 + - 53 ) COsb+ 2 ’, 5i + 52 < « < 

[ U - (gj + 52 + 53 ) + 5i + 52 + 53 < U < 



Then the cutter surface Ei is given by 


1 



MllM] 



Zcf 



Figure 3*1. Definition of a Generic Cutter 
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3.2.3 ALGORITHM: 

In thi.s algorithm denote the x, y coordinate values of the point J. 

1. Input the values of is-fillet-radius, d, r, a, b, hy 

2. If(is-fillet-radius){ 


d (cosa — sinb) 

e = r( -) 

2 cos{a + b) 

^ r + e sina 
cosa 

} 

else input the values of e and f 

3. Calculate the point P 


rP - 


= (e + / tana)cos^a — \/[(e + / tana)cos'^a]^ — (e^ + P — r'^)cos'^a 


4. Calculate the point Q 


= xf tana 


zf = {/ + [e + -{tana tanb — l)]tanb}cos^b 


+ 




[{/ + [e + ^{tana tanb — l)]fan6}cos^6]^ 
— {P -h [e + ^(tana tanb — 1)]^ — r'^]cos‘^b 


= zf tanb — 


-{tana tanb — 1) 


5. Calculate h 

If(is-fillet-radius) and (b = 0) h = zf 
else input the value of h 
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6. Calculate the point S 


Slid id 

fj. = n tanb — - iana tanb + - 

Jmt 

zf = h 


7. Calculate ^i,«S2,53 and $4 


= \/(a:fP + (zf)^ 


52 = r {arctan( 


/ 


Q 

Xc — e 


) — arctan(- 


/ 


)} 


.3 = \/(^f - *?)’ + (Zf - z?)" 


54 — 


8. Define the parametric equation '?(n) of the two dimensional profile of the 
cutter in terms of the parameter u 


, (it) = [a:c(it) yc{u) Zc{u) 1] 


as follows: 

yc{u) = 0,5i < u < 54 

lf(0 < u < Si){ 

Xc(it) = u cosa 
Zc{u) = u sina 


} else if(5i < u < $1 + S2){ 


. s . . ^ .sisina — f. ,ii — Si,. 

Xc(u) = e + lr| a> 5 {arctan(-^ — ^77;; p) + ( — — )} 


2c(w) = / + kl 5in{arctan( 


5i cosa — e r 

Si sina — f . — 5i 


5i cosa — e ' 




} 
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else if( 5 i + .S2 < u < Si + S2 + S3){ 

Xc{xi) = u - (si -i- S 2 + S 3 ) sinb + a:’ 

Zc{u) = u - (si + 52 + S3) cosb + zl 

} 

else if(si + 52 + S 3 < IX < Si + 52 + 53 + S 4 ){ 

2-c(u) = I* 

Zc{u) = U — (si + S2 + S3) + z" 

} 

9. Define the sweep matrix QM] 

10. Calculate -(u,v) using 

3.3 METHODOLOGY: 

In the previous chapters and the previous sections, the algorithms for the con- 
jugate surfaces and the shape of the cutter are discussed. However, they involve 
symbolic computation and hence cannot be implemented directly. A symbolic 
computation program MACSYMA is used to define this task. This finds the 
condition of contact and the corresponding x, y and z values. This section de- 
scribes the algorithm used for generating the conjugate surfaces at each time 
frame, the trimming surfaces and the algorithm used to show the animation of 
conjugate surfaces. 

Algorithm 3.3.1 explains the in-betweening and the generation of the conjugate 
surfaces while algorithm 3.3.2 describes the generation of the trimming surfaces 
and finally algorithm 3.3.3 illustrates the way animation is modeled. 
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3.3.1 ALGORITHM; 

1. Input 24 parameters of movement a,j, &,j, Cij, 9ij, 4>iji ’’Pij hJ = 

2. Input the 9 parameters of the cutter is — fillet — radius^a,b,d, e,/, h,hi,r 

3. Calculate 61,52,53,54 using algorithm 3.2.3 

4. For time = 0 to MAXT { 
t = time/MAXT 

for i = 1,2 { 

U, — Ctil “h (^{2 ^il)t 

bi = 6,1 + (6i2 — bii)t 
Ci = Cii + (Cj2 — Cii)t 
Bi — Bi\ + (^t"2 ~ Bii'^t 

<f>i = fix + ifi2 — 4>il)t 

iPi = + (l/>,-2 - fii)t 

} 

} 

5. For each of the four segments of the cutter, namely. 

Si to Si+i ,i = 0 to 3 

calculate X 2 {u,v),y 2 {u,v),Z 2 {u,v),con{u,v) using the algorithm 2.3 

6. For i = 0 to 3{ 

for U = 5; to Si+i do { 

a ) Calculate the values of v which satisfy the condition con(u,v) = 0 
b ) For each v in (a) do { 

1 . calculate 

X = X 2 {u, v) 


y=y2{u,v) 
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2 = Zi{u, u) 

2 . Draw the surface with values x, y, z 

} 

} 

} 

The algorithm needs some some explanation. The in-betweening is done using 
linear interpolation of the 24 parameters of movement. Out of these, twelve 
parameters corresponding to the initial position and the remaining twelve cor- 
respond to the final position. These are linearly interpolated in time. This is 
ensured in step 4. In step 5, ar 2 (u,u), j/ 2 (u, u), 22 (u, u) and con(u,v) are calcu- 
lated symbolically using the algorithm 2.3 . Here symbolic computation pro- 
gram MACSYMA is used to generate X 2 (u,v),y 2 (u,v),Z 2 (u,v) and con(u,v) in 
symbolic form. 

In reality, either all values of v satisfy the condition of contact con(u,v) in the 
range [0, 27 r] or there will be exactly two values of v satisfying the condition of 
contact in the range [0,27r]. Moreover, these two values of v will be almost tt 
apart. This is because the cutter shape is assumed to be axisymmetric. These 
values of v can be determined by a numerical method. ( In this thesis regula-falsi 
method is used). These give rise to eight different surfaces. 

Once the surfaces are drawn, the next step would be to trim the ends of the 
surfaces. This can be accomplished by drawing two more surfaces at the initial 
and final positions. The following algorithm illustrates the methodology to 
draw the trimming surfaces. An instance of the initial trimming surface and an 
instance of final trimming surface are shown in Fig.3-3. 

3.3.2 ALGORITHM: 

1. Put t = 0 in algorithm 3.3.1 

2. For each u in step 6 of the algorithm 3.3.1 do 

3. If there are not more than two roots in [0, 27r] for con(u,v) = 0 then 
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a) set Vi = root nearer to 0 

b) set V 2 = root nearer to vr 

4. For V = t?i to V 2 do { 

a) Calculate x, y and z using algorithm 3.2.3 

b) set X = X cost! 
set z = z sinv 

5. Plot X, y, z 

} 

} 

} 

6. Put t = MAXT in step 1 and repeat the steps 2 to 5 with v = t ;2 to vi •+ 27r 
in step 4. 

The last step is to get the trimming surface at the final position of the cutter. 
The first five steps give the initial trimming surface. It is seen that there will 
be either exactly two values of v or all values of v satisfying the condition of 
contact in the range [0,27r]. The trimming is not necessary when all the values 
of V satisfy the condition of contact. In other words, trimming needs to be done 
when there are only two values of v satisfying the condition of contact. Let the 
values be Ui and V 2 , v\ be the root nearer to 0 and V 2 be the root nearer to tt. 
The cutter shape being axisymmetric, the points need to be generated for the 
surface are nothing but the circular sweep from Vi to U 2 - The initial trimming 
surface is obtained when t = 0 and the final trimming surface is obtained when 
t = MAXT with the circular sweep being taken from V 2 to Uj + 27r. 

These algorithms provide the basis for drawing the surfaces. Now the following 
algorithm illustrates the displaying of the surfaces along with the cutter, in 
general, for each time frame. Some instances of the animation are shown in 
Figures 3-4 and 3-5. 
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3.3.3 ALGORITHM: 

1 . Draw the initial trimming surface using the algorithm 3.3.2 

2. For t = 0 to MAXT do { 

a) Draw the cutter surface using algorithm 3.2.3 

b) Draw the conjugate surfaces using algorithm 3.3.1 

} 

3. Draw' the final trimming surface using algorithm 3.3.1 

These algorithms are implemented in C on HP-9000/834 TURBO SRX 3D 
graphics workstations based on HPPA RISC CPU rated at 2.02 MFLOPS (14 
MIPS) performance with 24 plane 19” color monitor using the starbase graphics 
library. Smooth and continuous motion is achieved using the technique of double 
buffering in starbase. 












Chapter 4 


SOFTWARE DEVELOPMENT 


All the algorithms presented in this thesis are implemented on HP-900/834 
Turbo SRX 3D graphics workstation. The total package is divided into six 
modules (see Fig.4-1). Fig.4-1 illustrates the total methodology adopted for 
software development. 

4.1 CUTTER PROFILE MODULE : 

This module takes the nine parameters, namely, is-fillet-radius, a, b, d, e, f, 
h, hi and r. The module generates the two dimensional profile of the cutter 
according to the input parameters and develops the required input for the next 
module, namely, surface generation module. 

4.2 SYMBOLIC COMPUTATION MODULE: 

The symbolic algorithm is implemented using symbolic computation software 
MACSYMA (see Appendix A). This module requires the twenty four parameters 
in symbolic form. It generates the expressions for y 2 iu,v,t), Z 2 {u,v,t) 

and the condition of contact in terms of u, v and t. The same MACSYMA soft- 
ware may be used to eliminate one of u and v using condition of contact to 
get the surface Ej in terms of t and one of the parameters u and v. However, 
there are some inherent properties of the conjugate surfaces which are useful in 
implementation. For example, for each particular instance of u and t, either all 
values of v or exactly two values of v will satisfy the condition of contact in the 
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interval [0, 27r]. Moreover, if there are only two values of v satisfying the condi- 
tion of contact, then these two values of v differ almost by tt. This information 
is highly valuable in speeding up the symbolic program. This can be achieved 
by first identifying whether there are exactly two roots and then these two roots 
can be calculated faster using a numerical method. In this module the method 
of false position is used to calculate the values of v. 

The cutter 2D profile is axisymmetric and it consists of six straight line seg- 
ments and the two circular segments. Hence, the conjugate surface T>2 will be a 
combination of eight surface segments. Thus, it is sufficient to find the expres- 
sions for y2{u,v,t), Z2{u^v,t) and the condition of contact for any 

generic straight line segment and any generic circular segment. 

This module produces the expressions for X2(u,u,t), y2{u,v,t), Z2{u,v,t) and 
the condition of contact for any line segment with end points at (/i,?ni) and 
(/2,m2). Further, it produces the expressions for i2(w,v,t), y2{u,v,t), Z2{u,v,t) 
and the condition of contact for the circular segment, given the center and 
radius. These will be input to the surface generation module. 

4.3 SURFACE GENERATION MODULE: 

This module takes the twenty four parameters of motion and also takes the 
expressions for Z2(u, u, t), 7/2(11, u,t), 22(71, u,t) and the condition of contact 
in symbolic form. This module calculates the values of v using a numerical 
method ( the method of false position ) to eliminate v using the condition 
of contact. After v is eliminated from X2{u,v,t), y2{u,v^t), Z2{u,v,t) using 
condition of contact, the points X25 y2, ^2 are calculated. These Z2, y2i ^2 
account for the eight surfaces which form the conjugate surface S2. These data 
files ( generated by this module ) will be input for the trimming algorithm a^i 
well as the animation module. In fact, this module is the implementation of the 
algorithm 3.3.1. 
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4.4 (’UTTER SURFACE GENERATION MODULE: 

The cutter surface is a circular sweep of the 2D profile of the cutter that is 
(k'veloped by the cutter profile module. Hence, this module takes input from 
th<‘ cutter profile module and generates the cutter surface by taking circular 
sweep. Moreover, this module displays the cutter surface as a wire frame model 
and as a solid model. 





Figure 4- 1 Software Development 
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•1.5 THIMMINC; surface: generation MODULE: 

'Friimniiig is the* technique of generating the portion of the cutter surface that 
will he conjugate to the surface E 2 at the initial position of the cutter at t = 0 
and at the final position of the cutter at t = MAXT. 

This module takes the input from the surface generation module. The input will 
be the two values of v which satisfy the condition of contact at t = 0 and t = 
MAXT. Suppose iq and are the values of v satisfying the condition of contact 
at t = 0 for a part icular value of u. A semi-circular sweep from uj to V 2 or uj to 
I'l { depending on the orientation ) for all the values of u, gives the trimming 
surface* at the initial position. The trimming surface at the final position of the 
cutter can be generated in the same way except that the semi-circular sweep is 
to be taken frorxr rj to t'l ( or Ui to Uj depending on whether the initial trimming 
surface is taken from vi to V 2 or V 2 to wi). This module generates the trimming 
surfaces at the initial and the final positions of the cutter. This module is the 
implementation of the trimming algorithm 3.3.1. 

4.6 ANIMATION MODULE: 

'Fhe animation module is the implementation of the algorithm 3.3.3. This mod- 
ul<* draws the cutter surface at the initial position and with each time frame the 
cutter surface draws the corresponding conjugate surface and finally stops at 
the final position. The trimming surfaces are also drawn along with the cutter 

surface. 

This module is implemented using the double-buffering technique for movement 
and the light-source and rendering primitives of starbase graphics library for the 
rendering and shading the surfaces. 



Chapter 5 


CONCLUSIONS 


The work presented in this thesis outlines a methodology for displaying the con- 
jugate surface generation in animation. A generic cutter with all the twenty four 
parameters of movement is taken to illustrate the model. Linear interpolation 
has been used to give the effect of in-betweening in animation. 

The symbolic algorithm is used repeatedly to generate the conjugate surfaces 
and also for the animation purposes, Though, the work shown here is aimed at 
animating the conjugate surface generation, the algorithms and the methodol- 
ogy proposed in the thesis are general and can be used for modeling any shape 
gtmerating process. 

SCOPE FOR FUTURE WORK: 

The in-betweening technique used in this thesis is the linear interpolation. How- 
ever, in general, various interpolation schemes may be used depending on the 
application [79]. The interpolation schemes like spline interpolation are used 
in some applications of animation [51] and they may be used to enhance the 
speed and reality of this work. The method of quaternions proposed in [33] 
may improve the algorithms presented in this thesis, which might turn out to 
be faster than the conventional matrix multiplication methods. Most of the 
graphics softw'are is developed keeping in view the rendering and other related 
aspects. However, most of these are not suitable to applications of animation 
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systems, 'riiis may be due to tlie high data requirements of any animation sys- 
tems. Though there are quite few animation Ianguages[57,61,94], the problem 
needs to he addre.ssed in much more detail. 

'riu'oretically. if one can specify the class of all surfaces that can be generated 
using this model then it may be possible for the design engineers to choose the 
corresponding cutters and shapes. Another important aspect in the theoretical 
stand point is the limitations of the system and its practical implications. 




Appendix A 







MACSYMA 


Computers have traditionally been used to solve scientific problems that could 
be expressed in terms of numbers. When a problem can easily be expressed 
in terms of calculations with numbers, this approach to problem solving works 
well. On the other hand some problerns can be expressed best in symbolic 
terms or perhaps can only be expressed that way. MACSYMA is a symbolic 
algebra system. It can work with symbols, polynomial expressions, equations 
and numbers. 

M ACS^'MA is an interactive symbolic algebra program that offers a wide range 
of capabilities for solving many kinds of problems in algebra, trigonometry and 
calculus. The capabilities of MACSYMA include differentiating and integrating 
expressions, solving equations, manipulating matrices etc. MACSYMA is also 
a programming environment in which one can define mathematical procedures 
tailored to his own needs. MACSYMA can return results either in numeric or 
symbolic form. 

CAPABILITIES : 

a) The precision in the numerical systems is limited by the systems hardware. 
MACSYMA can work with exact quantities rather than approximations by 
carring out the computations in symbolic form. When a symbolic result is 
converted to a floating-point number, the precision can be set by the user. 


51 


b) (’oiislants like it and <• can be kept in symbolic form and use them even in 
floating-point calculations to get exact calculations using them. 

c) MAC'SYMA provides six categories of numbers, namely, INTEGERS, RA- 
TIONAL NUMBERS, FLOAT and BIGFLOAT and COMPLEX NUMBERS. 
MACS\'M.A does not limit the number of digits in an integer or a rational 
number but there are some limitations on the non-zero floating point numbers. 

(1) MACS^’MA can be viewed as a programming language. It accommodates 
many programming structures including conditional statements and loops. 

e) MACS^’MA provides the facility to convert any expression in symbolic form 
to the corresponding expression in C or FORTRAN or PASCAL. 



Appendix B 


STARBASE GRAPHICS LIBRARY 


Starba 3 e graphics library is a start of the art package of graphics procedures 
which, while being based on and supporting graphics standards, provides consid- 
erable functionality beyond that defined in graphics standards. The procedures 
defined in starbase can be accessed from the C, FORTRAN?? and PASCAL 
programming languages. Starbase is available on many HP-UX systems, in- 
cluding HP 9000 series 300 and 800 computers. Starbase is device independent. 
This appendix gives a very brief idea of the starbase graphics primitives. 

DOUBLE BUFFERING: 

Double buffering is a technique for increasing the smoothness of the drawing 
process on the screen. Double buffering involves taking a multiple-plane frame 
buffer ( 24 planes on TSRX systems ) and dividing it into two equal parts, Us- 
ing the double-buffer and dbufferjswitch procedures, one can ensure that 
only completed images are displayed. Once an image is complete, it is then 
displayed. The image in the other half of the frame buffer can then be mod- 
ified. Toggling back and forth in this manner is called double buffering. Any 
movement displayed using double buffering will be smooth and avoids jerkyness 
in movement. 

TRANSFORMATION MATRIX STACK: 

A transformation matrix performs scaling, rotation, translation etc on the co- 
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ordinate data. Starbase provides transformation matrices that one can control 
directly to suit the application . These matrices are arranged in a stack called 
transformation matrix stack. There are three types of transformations that can 
modify the stack. They are 

a) The modeling transformations 

b) The viewing transformations 

c) The VDC-to-DC transformations 

The modeling transformations are generally useful for the pre- multiplication /post- 
multiplication of the transformation matrices, pushing on to the stack, poping 
from the stack etc. some of the generally used routines of modeling transfor- 
mations are concat.transformationSd, pop_matrix, push_matrix3d etc. 

The viewing transformations control the synthetic camera, view-port, view- 
volume etc. The routines useful are view -Camera, view _port, view .volume, 
view .matrix . 

The VDC-to-DC transformations are useful in transforming from virtual de- 
vice coordinates to device coordinates. The routines useful in these are map- 
pingjriode, set. pi.p 2 ,vdc_extent, vdc.justification. 

HIDDEN SURFACE REMOVAL : 

The starbase uses z-buffer algorithm for hidden surface removal and can be 
accessed through the primitives hidden jsurface and zbuffer .switch. 

RENDERING: 

Starbase provides primitives for drawing wire frame models, shading and render- 
ing of the objects. Two primitives quadrilateral-mesh and triangular-mesh 
are useful for wire frame models. The light sources and illumination of am- 
bient, specular and reflections can also be modeled using the shade-mode, 
light -source and the surface-coefRcients routines. Starbase provides for 
both Phong and Garaud shadings. 
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Starliasc also provides primitives for modeling objects using B-spline curves and 
surfaces. Some more additional features like panning and zooming, displaying 
transparency and highlighting are also available in starbase. 
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